Infinite family of persistence exponents for interface fluctuations 
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We show experimentally and theoretically that the persistence of large deviations in equilibrium 
step fluctuations is characterized by an infinite family of independent exponents. These exponents 
are obtained by carefully analyzing dynamical experimental images of Al/Si(lll) and Ag(lll) equi- 
librium steps fluctuating at high (970K) and low (320K) temperatures respectively, and by quanti- 
tatively interpreting our observations on the basis of the corresponding coarse-grained discrete and 
continuum theoretical models for thermal surface step fluctuations under attachment /detachment 
("high-temperature") and edge-diffusion limited kinetics ("low-temperature") respectively. 
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PACS numbers: 68.35.Ja, 05.20.-y, 0.5.40.-a, 68.37.Ef 



Persistence is a fundamental and powerful concept in 
the stochastic dynamics of non-Markovian statistical pro- 
cesses Recently, the persistence concept has been 
applied to statistical studies of spatially extended non- 
equilibrium systems both theoretically and experi- 
mentally [3j. Loosely speaking, persistence is the prob- 
ability P(t) that a stochastic variable, which is found 
to have a "positive" value at time to, will stay "posi- 
tive" throughout a time interval t up to time to + 1. In 
most steady-state ("stationary") situations, the persis- 
tence probability is found to decay in time in a power- 
law manner, P(t) ~ t~ e , for large t, where the persis- 
tence exponent 8 is a highly nontrivial exponent charac- 
terizing the relevant stochastic dynamics under consider- 
ation. The concept of persistence has been used jiUHE] 
to study the interesting problem of thermally fluctuating 
interfaces where steps on vicinal surfaces undergo ran- 
dom thermal motion in equilibrium [7J. The step per- 
sistence probability is the probability P(to,to + t) that 
a given lateral step position x with a height (i.e. step 
fluctuation measured from the equilibrium step position) 
h(x,to) at time to does not return to this value up to a 
later time to + 1. With no loss of generality we will set 
to = from now on, assuming that thermal equilibrium 
has been achieved in the step fluctuations and we are dis- 
cussing steady-state stationary properties. The result- 
ing persistence probability P(t) = P(0, t) has recently 
been studied experimentally using dynamical scanning 
tunneling microscopy (STM) in two systems: Al steps 
on Si (111) surface at high temperatures ~ 970K Q and 
Ag (111) surface at low temperatures ~ 320K |fg. In 
the first case, the step fluctuations dominated by atom- 
istic attachment and detachment (AD) at the step edge 
are known |jj to be well described by the coarse-grained 
second-order non-conserved linear Langevin equation 



dh(x,t) _ d 2 h{x,t) 



+ j](x,t), 



(1) 



is the usual uncorrelated random Gaussian noise cor- 
responding to the non-conserved white noise associated 
with the random AD process. Low temperature step 
fluctuations dominated by the step edge diffusion (ED) 
mechanism are, on the other hand, described by a fourth 
order conserved linear Langevin equation: 



dh(x,t) _ d 4 h(x,t) 



dt 



dx 4 



Vc(x,t), 



(2) 



dt dx 2 
where rj[x,t) with (rj(x,t)r](x , t )) oc S(x — x )S(t — t ) 



where r) c (x,t) with (r] c (x,t)r) c (x , t )) oc V 2 5(x — x )S(t — 
t ) is a conserved noise associated with atomic diffusion 
along the step edge. From a quantitative analysis of the 
digitized STM step images as a function of time, the 
persistence exponent was found to be 8 — 0.77 ± 0.03 
P and 8 = 0.87 ± 0.02 @, respectively, for the high- 
temperature AD mechanism and low-temperature ED 
mechanism. These measured step fluctuation persistence 
exponents agree reasonably well with those found 4] from 
kinetic Monte Carlo simulations of the corresponding dis- 
crete solid-on-solid models: 8 ~ 0.75 (for Eq. JTJ) and 
8 ~ 0.86 (for Eq. (J5J). These results are in agreement 
with a postulated (and numerically verified) relation 0] 
between persistence in the steady state and dynamic scal- 
ing H in the pre-stationary transient regime. It is be- 
lieved that, at least for linear Langevin equations, the 
persistence exponent 8 is equal to (1 — (3) where (3 is the 
exponent that describes the initial power-law growth of 
the interface width Q in the transient regime ((3 = 1/4 
for Eq. and = 1/8 for Eq. ©). 

This striking agreement between experimentally ob- 
tained and theoretically predicted values of the persis- 
tence exponent demonstrates the overall excellent con- 
sistency among theory, experiment, and simulations in 
this problem, but also brings up a key question regarding 
persistence studies 0,0,0 of surface fluctuations: Is per- 
sistence really an independent (and new) conceptual tool 
in studying surface fluctuations, or, is it just an equiva- 
lent (perhaps even complementary) way of studying dy- 
namic scaling 0, of height correlations? In this Letter 
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we present new theoretical and experimental persistence 
results on surface step fluctuations that fundamentally 
transcend any dynamic scaling considerations, establish- 
ing in the process the existence of a novel and nontriv- 
ial infinite family (i.e. a continuous set) of persistence 
exponents for equilibrium step fluctuations. We carry 
out quantitative analyses of (digitized) dynamical STM 
images of step fluctuations both for high-temperature 
(Al/Si) and low-temperature (Ag) equilibrium surfaces, 
and compare in details the experimental results with 
those we have obtained from numerical integration of the 
corresponding Langevin equations and discrete stochas- 
tic Monte Carlo simulations of corresponding atomistic 
cellular automata type models in the same dynamical 
universality classes jjj 0. All three sets of persistence 
results agree very well for both high and low temper- 
ature equilibrium step fluctuations, establishing persis- 
tence (particularly, the infinite family of persistence ex- 
ponents) as a potentially powerful tool (rivaling, perhaps 
even exceeding, in utility the well-studied dynamic scal- 
ing approach) in studying dynamical interface fluctuation 
processes. 

The infinite family of persistence exponents we study 
here is based on the concept of persistence of large devia- 
tions introduced recently by Dornic and Godreche |3| in 
the context of kinetic Glauber-Ising dynamics of magneti- 
zation coarsening (a closely related idea, that of sign-time 
distribution, was developed in Ref. |l0|). For equilibrium 
step fluctuations, we define the probability of persistent 
large deviations, P(t, s), as the probability for the "aver- 
age sign" S'av of the height fluctuation to remain above 
a certain pre-assigned value "s" up to time t: 

P{t, s) = Prob { S av {t') > s, W < t }, (3) 

where 

S av (t) = t- 1 [ dt' Sit'), (4) 
Jo 

and 

S(t) = sign [h{x, t + 1) - h(x, t )], (5) 

where with no loss of generality, we set t = 0. Al- 
though the dynamical variable S(t) above is defined 
for a particular lateral position x, we take a statisti- 
cal ensemble average over all lateral positions to ob- 
tain a purely time dependent dynamical quantity P(t, s). 
Since S av (t) e [—1,1], the probability P(t,s) is defined 
for —1 < s < 1. For s — 1 we recover our earlier 
simple definition of persistence used in Refs. 0, 0, @ : 
P(t) = P(t, s = 1) measures the probability of the height 
fluctuation remaining above zero ("positive") through- 
out the whole time interval. For s = — 1 the probability 
P(£, s — — 1) is trivially equal to unity for all t. The gen- 
eralized probability function, P(t, s), defined in Eqs. J3J 
|J5J above, leads to a continuous family or hierarchy of 



persistent large deviations exponents, 0i(s), provided the 
steady-state decay of P(t, s) in time follows a power law, 
P(t, s) ~ t~ e '( s \ As we show below, this indeed happens 
for equilibrium step fluctuation phenomena, allowing us 
to define and measure the non-trivial persistence expo- 
nent 6i(s), —1 < s < 1, that varies continuously between 
8i(s = —1) = and 9i(s = +1) = 9 where 9 is the usual 
persistence exponent. Clearly, P(t, s) and 9i(s) are natu- 
ral generalizations of the persistence probability P(t) and 
the persistence exponent 9, respectively, to the broader 
concept of distribution of residence times with limiting 
behavior (i.e. s = 1) determining the usual persistence 
exponent. 

In this Letter, we have carried out the first application 
of the persistent large deviations concept to the equi- 
librium step fluctuations phenomenon. We also report 
the first experimental measurements of P(t, s) and 9i for 
any stochastic system. Our results for the two distinct 
types of step fluctuation processes (the high-temperature 
AD and low- temperature ED situations, described Q by 
Eqs. JU and J5J), respectively) are shown in Figs. d and 
The details of the experimental procedure for extracting 
the persistence probability from the digitized dynamical 
STM images are described in Refs. [a E| f° r the same 
two systems, and are therefore not repeated here. We 
only mention that in order to obtain the probability of 
persistent large deviations, we had to use a substantially 
larger number (~ 100) of STM data sets than that used in 
Refs. 5, 6]. As mentioned earlier, we compare the exper- 
imental results with two theoretical models in each case: 
The continuous Langevin equation (Eq. (JJJ or J2J) and 
a discrete stochastic model which is theoretically known 
[! El 0] to belong to the same universality class (in 
the dynamic scaling sense) as the continuum equation. 
For the high-temperature AD step fluctuations, described 
on a coarse-grained scale by Eq. Q (sometimes referred 
to as the Edwards- Wilkinson equation the discrete 
stochastic model we use is the extensively studied solid- 
on-solid Family model 0, 0] . For the low-temperature 
ED case, described on a coarse-grained scale by Eq. 
(sometimes referred to as the Mullins-Herring equation 
with conserved noise ||), the discrete stochastic model 
we use is the well-studied Racz model Our results 
are obtained from numerical integration of the Langevin 
equations using the simple Euler scheme Q, and stan- 
dard Monte Carlo simulation of the atomistic mod- 
els. Typical system sizes used in the numerical work are 
~ 1000 for Eq.JQ and — 100 for Eq.(J2J), and the number 
of independent runs used in the calculation of persistence 
probabilities is ~ 1000. 

In Fig. ^ we show in the top two panels our measured 
(panel (b)) and calculated (panel (a), main figure and 
inset) persistent large deviation probability P(t, s) as a 
function of time t for the high-temperature step fluctua- 
tions case. Each panel shows eight different log-log plots 
of P(i, s) against t for eight different values of the average 
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FIG. 1: Log-log plots of the persistent large deviations 
probability, P(t,s), for high-temperature surface step fluc- 
tuations via the AD mechanism, shown as a function of time 
(t) for different values of the average sign constraint param- 
eter s = 1.0, 0.75, 0.5, 0.25, 0,-0.25,-0.5,-0.75 (from the 
bottom to the top), (a) The continuum Langevin equation of 
Eq. Q (main figure) and the stochastic discrete Family model 
(inset); (b) experimental data from digitized dynamical STM 
step images of Al on Si(lll) surface; and (c) comparison of the 
various sets of results for 9i as a function of s. The error bars 
shown for the experimental data are obtained from variations 
of the local slope of the log P(t,s) vs \ogt plots. Simulation 
results for two sample sizes (L — 100 and L — 1000) are 
shown to illustrate that the use of small samples leads to an 
underestimation of 9i(s). 



sign parameter s (= +1, +0.75, —0.75 from the bot- 
tom to the top). As mentioned before, the case s = +1 
(the bottom-most curve) corresponds to the usual persis- 
tence probability P(t), and therefore the results shown in 
Fig.QJfor s = 1 are already known. The results for all the 
other values of s are new and non-trivial. The linearity 
of the log-log P(t, s) vs. t plots immediately implies that 
P(t,s) ~ t^ 6 '^. The different sets of results for 9i as a 
function of s are shown in Fig.^c). The excellent agree- 
ment among the various data sets shown in Fig. ^c) is 
the most important new quantitative result of our work. 
This means that the high-temperature step fluctuation 
phenomenon via the AD mechanism is indeed described 
by the Edwards- Wilkinson equation (and therefore also 
by the discrete Family model), not just in the sense of the 
dynamical universality class (as defined by specific expo- 
nent values, e.g. f3 and 9) but more importantly for the 
infinite family of persistent large deviations exponents as 
defined by the continuous function 9i(s). This striking 
agreement between experiment and theory for a contin- 
uous family of exponents definitely establishes persistent 
large deviations studies as a new and effective tool for 
studying dynamical fluctuations of nanoscale systems. 

In Fig. |2 we present results similar to those in Fig. ^ 
but now for the ED mechanism step fluctuation data 
along with the corresponding theoretical results for the 
continuum Langevin equation defined by Eq. (J2J) and 
the discrete stochastic Racz model which are known 
0, IE lul t° be in the same dynamic universality 
class as the low-temperature step fluctuation process. 
The same description and explanation given above for 
Fig. □ apply now to Fig. where 9i(s = 1) « 0.875, 
which agrees with recent experimental measurements 
of the usual persistence exponent of low-temperature step 
fluctuations on Ag and Pb surfaces. The experimental 
and theoretical results for the continuous function 9i(s), 
shown in Fig. Efc), exhibit qualitative agreement, with 
the experimental exponent values for s < 1 being slightly 
larger than the theoretical ones. There are several pos- 
sible explanations for this difference between experimen- 
tal and theoretical results. There are reasons to expect 
that increasing the dynamic range of the experimental 
P(t, s) beyond two decades in t (this range is limited 
by noise problems inherent in dynamic STM imaging) 
would bring theory and experiment into closer agree- 
ment. To illustrate this possibility, we show in the inset 
of Fig. |2{c) the time-dependence of the local exponent 
9i(s,t) = d \og[P (s,t)]/d log t, obtained from simulations 
of the Langevin equation for two values of s for which 
the difference between theory and experiment is large. 
The local exponent is found to decrease with time before 
reaching a constant value at large t. We have checked 
that the experimental data show similar behavior for all 
s < 1, which implies that the effective exponent values 
obtained from power-law fits of relatively short-time data 
would be larger than the true long-time values. Indeed, 
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FIG. 2: The same as in Fig. Q for the low-temperature 
step fluctuations via the ED mechanism, (a) The contin- 
uum Langevin equation of Eq. @ (main figure) and the dis- 
crete stochastic Racz model (inset); (b) experimental data 
for the Ag(lll) step fluctuations; and (c) comparison of the 
three sets of results for 9i : experiment (star) , stochastic Racz 
model (square) and Eq. JS) (circle). The inset shows the time- 
dependence of the local exponent 8i(s,t) obtained from simu- 
lations of the Langevin equation for s = 0.25 and s — 0.50. 

we have found that fits of the simulation data over the 
range 20 < t < 200 (this is the range used in obtaining 
0i(s) from the experimental data) yield values of 0i(s) 
that are higher and closer to the experimental values. 
A second possibility is that the smallness of the sample 
size used in the simulations leads to an underestimation 
of the values of 6i(s), as in Fig. ^c). Unfortunately, 
the impossibility of equilibrating much larger samples of 



these models with very slow dynamics prevents us from 
checking this explicitly. Another possibility that we can 
not rule out is that the Ag(lll) equilibrium step fluctu- 
ations do not precisely follow the theoretical models of 
edge-diffusion limited kinetics. Further experimental and 
theoretical investigations would be needed for settling 
this issue. We should emphasize, however, that given 
the severe complexity in measuring any power law expo- 
nents associated with surface step fluctuation dynamics, 
the overall agreement between theory and experiment is 
quite good. 

In summary, we have established the concept (and 
the usefulness) of an infinite family of persistent large 
deviations exponents for height fluctuations in equilib- 
rium surface step dynamics phenomenon. The impres- 
sive agreement between theory and experiment indicates 
that the persistent large deviations probability (and the 
corresponding exponents) may very well be an extremely 
powerful tool in characterizing and understanding other 
stochastic fluctuation phenomena, e.g., kinetic surface 
roughening in noncquilibrium growth. In contrast to 
other dynamical approaches, the technique developed in 
this Letter leads to a continuous family of exponents (a 
continuous function rather than one or two isolated inde- 
pendent exponents (as in the dynamic scaling approach) 
and is therefore a much more stringent test of theoretical 
ideas, and also perhaps provides a deeper level of probing 
the dynamics of fluctuation problems. 
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